Cold neutron scattering in imperfect deuterium crystals 

Andrzej AdamczaliE 

Institute of Nuclear Physics, Polish Academy of Sciences, 
Radzikowskiego 152, PL-31342 Krakow, Poland 
(Dated: December 13, 2010) 

Abstract 

The differential cross sections for cold neutron scattering in mosaic deuterium crystals have been cal- 
culated for various target temperatures. The theoretical results are compared with the recent experimental 
data for the neutron wavelengths A « 1-9 A. It is shown that the structures of observed Bragg peaks can 
be explained by the mosaic spread of about 3° and contributions from a limited number of crystal orienta- 
tions. Such a crystal structure should be also taken into account in ultracold neutron upscattering due to the 
coherent phonon annihilation in solid deuterium. 
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I. INTRODUCTION 



The aim of this work is a calculation of the differential cross sections for cold neutron scat- 
tering in solid molecular deuterium (sD 2 ) at various temperatures and comparison with the recent 
experimental data OLSl- Deceleration of cold neutrons in sD 2 converters is an effective method for 
production of high-intensity ultracold neutron sources OHsl]- Therefore knowledge of the neutron 
cross sections is indispensable for effective projecting such sources. Ultracold neutrons (UCN) 
are applied for studies of fundamental properties of the neutron [1(1 11]. Also, they can be used 
for surface physics 1120, investigations of the neutron mirror-neutron oscillations [13] and neutron 
quantum states in the gravitational potential II 1411 . 

The experimental total cross sections OL 0] in the Bragg-scattering region display a certain 
number of pronounced peaks. Their heights decrease when sD 2 temperature is rising from 8 K to 
18 K. However, for a fixed crystal, the structure of Bragg peaks does not practically change with 
varying temperature. It is well known that sD 2 at such temperatures and under low pressure has 
the hep structure [15, lq) . which has been confirmed using the Raman spectroscopy 0]. A simple 
calculation of the Bragg cross sections presented in Ref. LI] does not lead to agreement with the 
experiment, both for the random poly crystalline hep structure and for a few specifically oriented 
hep polycrystals. The magnitudes of the theoretical peaks are too large. Also, these peaks are 
sharp, while the experimental peaks are quite broad. Moreover, the observed and calculated lo- 
cations of the Bragg peaks do not fully coincide. A more advanced calculation of the total cross 
section for a random poly crystalline sD 2 II 1711 . with incoherent processes taken into account, also 
disagrees with the experiment in the Bragg-scattering region. 

A calculation of the cross sections presented in this paper includes both the elastic and in- 
elastic processes in the sD 2 lattice and in D 2 molecules. Rotational and vibrational transitions in 
these molecules, which are due to collisions with neutrons, are taken into account. Interactions 
of neutrons with phonons are described in the incoherent approximation while the incoherent and 
coherent elastic cross sections are treated separately. The coherent elastic scattering of neutrons is 
considered both for the random polycrystalline hep structure and for different orientations of the 
hep monocrystal. 

The calculated total cross sections are compared with the available data for the neutron wave- 
lengths A « 1-9 A |Q1 @], with the experimental uncertainty of neutron wavelength taken into 
account. Since this uncertainty is much smaller than the observed widths of the Bragg peaks [18], 
an explanation of this effect by the presence of mosaic imperfections in the sD 2 crystals is studied. 
Most crystals are perfect only in very small domains (mosaic blocks) on the order of 1000 A, 
which are separated by displacements and distortions in the lattice. These domains are almost par- 
allel, with small deviations of their orientations from the main direction. Such angular deviations 
are usually well-described by a Gaussian distribution 111911 . The mosaic spread may reach even 
a few degrees in certain crystals. 

Another experimental problem is connected with different structures of the sD 2 crystals, which 
are grown even in very similar conditions. As a result, the experimental Bragg cross sections for 
the sD 2 crystals, obtained by freezing D 2 gas and two samples of liquid D 2 , are different 01 Q]. In 
this paper, such a phenomenon is ascribed to different orientations of large monocrystals (much 
greater than the mosaic blocs) that are present in a given sD 2 target. These orientations differ much 
greater than those of microscopic mosaic blocs within a single monocrystal. Such a model of a 
real deuterium target is suggested by the optical observations of growing sD 2 01 @]. 
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II. METHOD OF CALCULATION 



Interaction of the neutron with a single deuterium nucleus is described using the Fermi pseu- 
dopotential [20]. The partial differential cross section for neutron scattering in sD 2 can be ex- 



pressed in terms of the Van Hove response function S [21], which depends on the target prop- 
erties for fixed momentum transfer and energy transfers. When the impinging neutron causes a 
rotational- vibrational transition in a deuterium molecule bound in the sD 2 lattice, the correspond- 
ing cross section can be expressed as follows 

d 2 a \ W 

= —O- 0n Si(K,CO), (1) 



\dQds') 0n k 



where <S, denotes the incoherent response function, which is a fraction of the total response func- 



tion S Il22h . The energy transfer co and the momentum transfer k to the sD 2 lattice are defined as 
(in the atomic units) 

to = s - s' - AE , k = k - k', (2) 

where s and s' stand for the initial and final kinetic energies of the neutron and AE is the internal- 
energy change of the target molecule due to the rotational- vibrational transition. Vectors k and k' 
denote the initial and final neutron momenta, respectively. The usual relations between the kinetic 
energies e and e' and the corresponding absolute values k and k! of neutron momenta are fulfilled 

e = k 2 /2m n , s = k' 2 /2m n , (3) 

m n being the neutron mass. Function cr „ in Eq. §\§ is the squared modulus of neutron scattering 
amplitude for a single bound molecule 

cr 0n = \r 0n \ 2 , (4) 

averaged over the total spin of the neutron and D 2 system. The indices (0) = (OK) and (ri) = (vK') 
denote the initial and final rotational- vibrational states, respectively. In the case of considered 
target temperatures, all the molecules are initially in the ground vibrational state. The final vibra- 
tional state is characterized by the quantum number v. The quantum numbers K and K' correspond 
to the initial and final rotational states of D 2 . 

When AE = 0, for all the molecules in sD 2 , neutron scattering may be strictly elastic (no 
simultaneous phonon processes) or quasielastic (with phonon creation or annihilation). It is con- 
ventional to separately calculate the coherent and incoherent fractions of the corresponding differ- 
ential cross section. The coherent cross section, which displays strong interference effects from 



all lattice sites, takes the general form II22I1 



d 2 cr \ k! 

— cr coh SiK, 0)), (5) 



per one D 2 molecule. The coherent elastic cross section cr coh for a single bound molecule is defined 

as 

C"coh = \Tqo I . (6) 

The horizontal bar denotes here averaging over a random distribution of the total spin of the n + D 2 
system and over the molecular rotational population in the lattice. The value / of total spin of the 
symmetric molecule D 2 is equal to / = 0, 2 for even rotational numbers K (ortho-D 2 ) and 1=1 
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for odd rotational numbers (para-D 2 ). In solid deuterium, the two rotational states K = and 
K = 1 are usually present, since the rotational relaxation K = 1 — » K' = is very slow without a 
catalyst 111 611 . It is assumed that there is no correlation between the spin I and the location of D 2 in 
the sD 2 lattice. 

The incoherent cross section for neutron scattering in sD 2 (A£ = 0) is expressed by the formula 



similar to Eq. © MM 



dl(T ^ k ' Q / \ m 
— o- inc Si(K,oj), (7) 



where cr inc is the incoherent elastic cross section for neutron scattering from a single bound 
molecule 

2 I 1 2 ^ 
C"inc = |!mX)I _ pOo| = I^Ool _ 0"coh • (8) 

The functions cr „, <x coh , and cr inc have been calculated similarly as the analogous quantities for 
the muonic atom scattering in gaseous and solid hydrogenic targets 11231 UM. Since the muonic 
atoms are small neutral objects, the methods of evaluation of the cross sections for these atoms 
are almost the same as for neutrons. Therefore, it is sufficient to use the neutron mass, scattering 



lengths, and appropriate spin correlations |25Ll260 in the general equations from Refs. [|23|,|24f]. As 



a result, in the case of n + D 2 scattering one obtains 



O"0n 



( 4 ^coh + 5 i b U fKK'(K) even K and K', 
§*L even K, odd K' , 

odd K, even K' , 



.( 4 ^coh + t>l c )f K K>(K) odd K and K', 

where = 6.671 fm and b[ nc = 4.04 fm denote, respectively, the coherent and incoherent 
scattering length for neutron scattering from a bound deuterium nucleus. The molecular form 
factor f KK > is defined as 

= YW K , IK \D VI (K)\ 2 , (10) 



with the angular- momentum factor 'Wk'Ik 



2 

<W mK = (2T + 1)(2/ + 1) ( o' o o) ' = 1 (1 1} 

expressed by the Wigner 3 j symbols. Equation (fTTT) has been obtained upon averaging <x 0n over 
projections of the initial angular momentum of D 2 and summing over projections of the final 
angular momentum of this molecule. Function T> v \ in Eq. (fTOl) denotes the radial matrix element 

dRx v (RM<<Rn)xo(R), 



(12) 



when v = and / = , 



*->o otherwise, 



in which j/ is the Zth spherical Bessel function. The radial wave function of the D 2 molecule 
for the vibrational state v is denoted here by^ v and R stands for the internuclear distance in D 2 . 
Numerical evaluation of the matrix elements (fT2l) has been performed assuming that the vibrations 
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of D 2 are harmonic and that there is no coupling between the molecular vibrations and rotations. 
Employing the harmonic -potential model of D 2 instead of the rigid-rotor model leads to accurate 
cross sections even at high collision energies. 

The coherent elastic cross section © is expressed by the coherent scattering length b coh and the 
radial matrix element (PT2l) for v = and / = 

K 2 oh l£>oo(*)l 2 for K = K', 
Ccoh = < n , . (13) 

10 otherwise. 

In the case of harmonic hep lattice, the incoherent response function Si in Eqs. CQ) and © is 
approximated by the phonon expansion of S t for a Bravais cubic lattice &2V Wh 



Si(K,cS) = exp(-2W) 



^ (2M 

8(lO) + ) gn(0j) 

n=\ 



(iwy 

! 



(14) 



in which exp(-2W) = exp[-2W(ic)] denotes the Debye-Waller factor. The functions g n are defined 

as 

1 K 2 Z(C0) 

2W 2M mo \ co 

gn(C0)= I dco'g^O)- 0)') gn-tico') , (15) 



-L 
f 



diOgnico) = 1 



where M mo i is the mass of D 2 molecule, Z(co) is the normalized density of vibrational states of the 
lattice, and n B denotes the Bose phonon-population factor 

n B (co) = [exp(a>/k B T) - l]" 1 , (16) 

for a given temperature T (k B is the Boltzmann constant). The first term in the phonon expansion 
(fT4l) corresponds to recoil-less scattering from a rigid sD 2 lattice. When AE = 0, such a scattering 
cannot change the neutron kinetic energy because of a very large mass of the sD 2 target. The 
next terms describe neutron scattering and possible rotational- vibrational transitions in D 2 with 
simultaneous creation (co > 0) or annihilation (co < 0) of phonons. In particular, the second term 
in Eq. (fT4l) is connected with one-phonon processes. Equations (TT5T) have been derived assuming 
thatZ(-w) = Z(co). 

Evaluation of the incoherent cross section has been performed using the isotropic Debye model, 
characterized by the following density of vibrational states 

Z M = M< V"* 00 (17) 

10 II CO > COb , 

where w D = k B D denotes the Debye energy corresponding to the Debye temperature D . In the 
case of solid hydrogenic targets at low pressures, 6>d equals about 100 K and it slowly decreases 
with rising target temperature [fl5ll . For the calculation of cross sections at various temperatures, 
the following estimation of the Debye temperature 

6> D = 2 1/2 (3tt 2 ) 1/3 /3 — ^— , (18) 

^hep 
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has been used. The sound velocity c s and the molar volume of sD 2 for calculation of the lattice 
constant ah cp have been taken from Ref. lll6ll . A constant correction factor fj = 1.085 in Eq. (PT81) 



is introduced to fit the experimental value of 1 14 K for sD 2 at zero temperature [28]. As a result, 
one obtains <9 D = 1 10 K for 8-K sD 2 and 6> D = 104 K for 18-K sD 2 . 



The exponent 2W of the Debye- Waller factor for a cubic lattice takes the form 112211 



2W=\{xx 2 )k 2 , (19) 



where (u 2 ) is the mean square displacement of the molecule from its lattice site. In the case of 
harmonic crystal, Eq. (fT9~l) can be estimated using the expression 

2W=-^- r^cothf-^)^. (20) 



2M mol Jo \2k B T ) co 

This expression is employed here as a good approximation for a harmonic hep lattice. For a 5-K 
sD 2 crystal, Eqs. (fT8l and ([201 ) lead to the root-mean- square displacement of 0.50 A, which is in 
excellent agreement with the experimental result from Ref. [28fl. 

The calculation of all phonon cross sections is performed here in the incoherent approach, 
which is achieved by substituting <x inc + cr coh for cr inc in the nonelastic fraction of Eq. ©. As a 
result, only the elastic Bragg scattering is taken into account in the coherent cross section ©. This 
is a very good approach at s > 1 meV. For a perfect lattice, the coherent elastic cross section per 
one molecule has the form of a sum over the reciprocal-lattice vectors t [22] 

f^) = X \ F N(r)\ 2 SCc - r) exp(-2W) , (21) 



where V denotes the unit-cell volume and \F n (t)\ is the average squared unit-cell factor 

N 

N I 



1 N 2 
F n (t)\ 2 = cr coh -|^exp(rr-d) . (22) 

d=i 

Vector d denotes here the position of a given molecule in the unit cell and N is the number of 
molecules per unit cell. The unit cell of hep sD 2 lattice contains two molecules. From Eq. (|2"TT) it 
follows that the Bragg scattering takes place only if the condition 

K = k - k' = r (23) 

is fulfilled. This condition cannot possibly be satisfied below the Bragg-cutoff energy 

s B = r 2 nin /(8m, 1 ), (24) 

in which T min denotes the value of the shortest vector t with F n (t) 0. In the case of hep crystal 
with the following basic vectors of the inverse lattice 

n =c(l,l/V3,0), r 2 = c(0,2/V3,0), 
t 3 = c(0,0, V3/8), 

where c = 2n/a hcp , the Bragg-cutoff energy is determined by the vector r m j n = [1,0,0]. This 
corresponds to e B ~ 2 meV for sD 2 at low pressure. 
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When SD2 has a random polycrystalline structure, the cross section (|2TT) . derived for a 
monocrystal sample, can be averaged over all orientations of the lattice. This leads to the poly- 
crystalline cross section 

§£" f ^l^^l-^-co S(? )exp(-2W), (26) 

# being the angle between vectors k and k'. If a mosaic sD 2 monocrystal is considered, Eq. (12TT) 
is numerically averaged over the Gaussian distribution 

g(<*) = — Wxp(-— )• (27) 



or 



V27T \ 2a, 



of angular deviation a of the mosaic blocks from the mean direction. The standard deviation aro, 
which is usually called the mosaic spread, is here a free parameter. 

The total differential cross section for neutron scattering in sD 2 is a sum of the incoherent cross 
sections (Q~|) and ©, with Si given by Eq. (fl4l) . and the coherent Bragg cross section(l2"TT). For 
fixed initial energy s and rotational state K, all possible rotational and vibrational transitions are 
included in this cross section. The total cross section cr tot is thus given as 



0~ to 



where the coherent fraction is calculated using Eq. ((21) or Eq. (1261) . depending on the structure 
of a sample. If several rotational states are initially populated, the incoherent fraction of cr tot is 
additionally averaged over a distribution of these states. At e » sb, the coherent effects are 
negligible and therefore the total cross section is solely determined by the incoherent processes. 
At s » Wd, effects of molecular binding in the lattice can be neglected and thus the cross section 
for sD 2 (per one molecule) corresponds to that for a free D 2 molecule. 



III. RESULTS 

Below, the calculated cross sections are compared with the experimental total cross sections 
for cold neutrons 01 0], which were determined from neutron-beam attenuation in sD 2 samples. 
The purpose of the experiment was to study sD 2 as the UCN converter. The measurements were 
performed in almost pure ortho-D 2 (98.7 + 0.2%), since a significant admixture of para-D 2 would 
drastically decrease the UCN lifetime rf29ll . due to the neutron energy gain during the rotational 
relaxation K = 1 — » K' = in D 2 . Therefore, the present calculations have also been performed for 
pure ortho-D 2 . In comparison with the experimental data it is assumed that the sD 2 target is small 
so that the neutron leaves the crystal after a single scattering and thus the kinematic diffraction 
theory is applicable. All available data are normalized to the value of 4.1 b at A = 6.5 A and 
T = 8 K [18], which is due to problems in accurate normalization of the measured cross sections. 

The calculated total cross section for cold neutron scattering in a polycrystalline sD 2 at tem- 
perature T = 8 K and 15 K is shown in FigQ] One can see that a very good agreement with 
experiment is reached for A > 6.5 A, where there is no Bragg scattering and the total cross sec- 
tion is practically determined by the incoherent processes. Also, a good agreement is visible for 
A > 2 A, where a contribution of the Bragg scattering to the total cross section is not dominant. 
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FIG. 1: Calculated (solid and dashed lines) and experimental (squares and circles, from Fig. 4 of Ref. 
total cross sections for cold neutron scattering in a polycrystalline ortho-D2 at temperature T - 8 K and 
18 K as functions of neutron wavelength. Also are shown the calculated incoherent contributions to the total 
cross sections (dotted and dot-dashed lines). 



The incoherent cross section is greater at T = 18 K since the phonon processes are more probable 
at higher temperatures. On the other hand, the calculated Bragg cross section is greater at T = 8 K 
because the Debye-Waller factor is smaller at higher temperatures. Similar results were obtained 
in Ref. y/ZD, with the use of Z(a>) derived from experiments. 

Different contributions to the total incoherent cross section are shown in Fig. [2l for a solid 
ortho-D 2 at T = 8 K. The interval A « 4-7 A is dominated by the incoherent strictly-elastic 
scattering. The cross section of this process is described by the first term of expansion (fl~4~l) . The 
purely phonon processes are most significant in the vicinity of A = 2 A. The rotational- vibrational 
excitations begin to dominate the incoherent cross section above the first rotational threshold at 
A « 1 A, which corresponds to the excitation K = — » K' = 1. Since in sD 2 this excitation 
may take place with simultaneous phonon annihilation, the respective cross section is visible up to 
A~ 3 A. 

The theoretical curves in Fig. Q] suggest that the experimental samples were not random sD 2 
polycrystals. In the interval A « 2-6.5 A, the experimental data display several broad peaks with 
similar widths. The magnitudes of theoretical sharp peaks are mostly much larger and practically 
do not depend on temperature. Their locations do not fully coincide with the locations of observed 
peaks. In order to understand this problem, it is reasonable to consider neutron scattering in sD 2 
monocrystal at various orientations. In comparison of theory with experiment, it is necessary to 
take into account the wavelength uncertainty AA « 0.05 A fl 1 8M . This value is too small to explain 
the presence of broad Bragg peaks in a monocrystal sample. Therefore, a possible mosaic structure 
in the sD 2 monocrystal is considered, which can lead to a significant broadening for a sufficiently 
large value of the mosaic spread. 

The calculated total cross section for cold neutron scattering from a 8-K sD 2 monocrystal ori- 
ented in the [2,0,1] direction is shown in Fig. |3] A mosaic structure of the monocrystal, with the 
spread chq = 3°, has been assumed. Also, the wavelength resolution is taken into account. The plot- 
ted experimental data [Ql 0] for the three 8-K samples demonstrate that the specific structures of 
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FIG. 2: Calculated contributions to the total incoherent cross section as functions of the neutron wavelength, 
for a solid ortho-D2 at T = 8 K. The label "elastic" corresponds to the incoherent strictly-elastic scatter- 
ing. A sum of the phonon creation and annihilation cross sections, with AE = 0, is denoted by the label 
"phonon". A sum of the rotational-vibrational transitions with possible simultaneous phonon processes is 
represented by the label "rovib". 



sD 2 were different, although freezing conditions were very similar. The two samples were frozen 
from liquid D 2 (labels "liquid 1" and "liquid 2") and one sample was obtained by freezing D 2 gas 
(label "gas"). The presence of various patterns of Bragg scattering is obvious, if the considered 
samples were, at least in part, monocrystals grown in different orientations with respect to the 
incident neutron beam. Some peaks of the calculated cross section approximately correspond to 
the peaks observed in the "liquid 2" target. On the other hand, the calculated Bragg maximum at 
4 A is visible in the "liquid 1" case. The mosaic spread of about 3° is optimal — an increase or 
decrease by 1° leads to a worse agreement with the experiment. As one may expect, in the limit 
ao — » 180° the theoretical cross section tends to the polycrystalline pattern, independently of a 
chosen crystal orientation. 

The cross sections calculated for all the orientations [h, k, I] with h, k, I = 0, 1, . . . , 10 do not 
accurately reveal the observed Bragg patterns of the three samples. This suggest that these samples 
could be composed of a certain number of large mosaic monocrystals. Visible cracks and optical 
opacity of the sD 2 targets [HI 01, which vary with temperature, support such a model. The total 
cross section in Fig.H]has been obtained for the superposition of two crystal orientations: [5,0,2] 
with the probability V = 0.45 and [5,5,2] with V = 0.55. The theory quite well reveals the 
locations of Bragg peaks for the "liquid 2" sample. The differences between magnitudes of the 
calculated and experimental peaks do not exceed 30%. In the "liquid 1" case, which is presented 
in Fig. |5l an approximate description of the data has been achieved assuming the presence of three 
following monocrystal orientations in the sample: [1,0,3], [1,1,3], and [2,0,1]. Here they have the 
same probability V = 0.33. 

The peaks of the measured cross section for the "gas" sample are less pronounced than the 
peaks observed in the "liquid 1" and "liquid 2" targets. This can probably be ascribed to the 
lack of local correlations between D 2 molecules in gaseous deuterium. As a result, much more 
various orientations of sD 2 monocrystals are present in the "gas" sample, compared to the targets 
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FIG. 3: Calculated total and incoherent cross section for cold neutron scattering from a 8-K ortho-D2 

monocrystal, oriented in the [2,0,1] direction, as a function of neutron wavelength. The mosaic spread 

n 

ao = 3° has been assumed. The experimental data for T = 8 K (from Fig. 5 of Ref. [2]) were determined 
using three different SD2 targets, obtained by freezing D2 gas (one sample) and D2 liquid (two samples). 




FIG. 4: Calculated total and incoherent cross section for cold neutron scattering from a 8-K ortho-D2 target, 
which is composed of large hep monocrystals oriented in the directions [5,0,2] (45%) and [5,5,2] (55%). 
The mosaic spread of 3° is assumed. The experimental 8-K data were measured using the SD2 sample grown 
from liquid D 2 ("liquid 2", Fig. 5 from Ref 



grown from liquid deuterium. The theoretical cross section in Fig.|6]has been calculated assuming 
that the "gas" target is composed of the monocrystals that have orientations corresponding to the 
69 shortest vectors of the inverse lattice (in the range [0,0,1]-[5,1,0]). It is assumed that all the 
orientations have the same probability. The mosaic spread of 3° in every single monocrystal is 
also taken into account. One can see that the calculated total cross section roughly corresponds to 
the "gas" data. A better agreement with the experiment can be achieved by tedious adjusting the 
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FIG. 5: Calculated total and incoherent cross section for cold neutron scattering from a 8-K ortho-D2 target, 
which is composed of large hep monocrystals oriented in the directions [1,0,3] (33%), [1,1,3] (33%), and 
[2,0,1] (33%). The mosaic spread of 3° is assumed. The experimental 8-K data were measured using the 
SD2 sample grown from liquid D2 ("liquid 1", Fig. 5 from Ref yfl). 




FIG. 6: Calculated total and incoherent cross section for cold neutron scattering from a 8-K ortho-D2 target, 
which is composed of large hep monocrystals oriented in all directions corresponding to the 69 shortest 
vectors of the inverse lattice (in the range [0,0,1]-[5,1,0]). The mosaic spread of 3° is assumed. The 
experimental 8-K data were measured using the SD2 sample grown from a gaseous D2 ("gas", Fig. 5 from 
Ref 101). 



probabilities of various orientations, which are free parameters, and by using more monocrystal 
orientations, especially in the "liquid 1" and "liquid 2" cases. However, this is not a purpose of 
this work. 

The experimental data shown in Fig. \T\ were measured using the same "liquid 2" sample at 
T = 8 K and 18 K. The Bragg patterns for different temperatures, including 14 K and 17 K [3,0], 
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are very similar. This means that the mosaic spread and orientations of macroscopic monocrystals 
do not significantly change, for a fixed sample. On the other hand, the magnitudes of observed 
Bragg peaks systematically decrease with rising temperature. This decrease is appreciable, es- 
pecially when compared to the incoherent cross section which increases with rising temperature. 
The incoherent processes are well described by the theory. However, the calculated Bragg cross 
sections do not reveal strong temperature changes of the measured cross sections, both for the 
random-polycrystalline and monocrystal cases. A possible explanation is that multiple Bragg 
scattering took place in the sD 2 targets and, thus, the kinematic diffraction theory is not fully 
applicable to the considered experiment. In such cases, phenomena of primary and secondary 
extinction, which strongly depend on crystal structure and change the neutron-beam attenuation, 
should be taken into account (see e.g., Ref. [30] and references therein). In Fig. [7} the calculated 




FIG. 7: The same as in Fig. 0]for T - 8 and 18 K of the "liquid 2" sample. A contribution from the Bragg 
scattering to the calculated total cross section at T - 18 K is scaled by the factor of 0.5. 



total cross sections for the "liquid 2" sample at 8 K and 18 K are compared with the experimental 
cross sections. A qualitative agreement of theory and experiment is achieved here by scaling the 
8-K Bragg cross section by a factor of 0.5. 



IV. CONCLUSIONS 



The calculated cross sections for cold neutron scattering in solid deuterium are in a reason- 
able agreement with the experimental data, if it is assumed that the sD 2 targets are composed of 
limited numbers of mosaic monocrystals with the mosaic spread of about 3°. In particular, this 
naturally explains the presence of different Bragg patterns in the sD 2 samples which are grown 
separately, although under very similar conditions. Such patterns are due to various orientations 
of the monocrystals, with respect to the impinging neutron beam. The magnitude of incoherent 
cross section and its variation with temperature is well described by the theory. On the other hand, 
the magnitude of measured coherent elastic cross section decreases strongly with rising tempera- 
ture, which is not revealed by the theory. A possible explanation of this disagreement is that the 
kinematic diffusion theory is not valid for comparison of the theoretical cross sections with the 
experimental cross sections, which are determined using the neutron-beam attenuation method. 
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Although the presented cross sections have been calculated in the Bragg-scattering region, the 
conclusions regarding the sD 2 structure are of importance for scattering of ultracold neutrons and 
thus for developing high-intensity UCN sources. The phonon annihilation in ortho-D 2 is a crucial 
process which leads to UCN upscattering and, therefore, limits the effectiveness of sD 2 converters. 
Since the amplitude of neutron scattering from deuterium is mostly coherent, it is necessary to 
take into account a real target structure in the UCN-energy region. In estimations of the coher- 
ent phonon annihilation, the polycrystalline approach and neglecting crystal imperfections may 
lead to strong disagreement between theory and experiment. Figure 2 of Ref. 113 ill , which shows 
strong changes of the UCN cross sections for different 5-K sD 2 samples, is an example of such 
disagreement. 
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